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Abstract 

We perform a comprehensive survey of the potential energy landscapes of 13-atom 
Morse clusters, and describe how they can be characterized and visualized. Our aim 
is to detail how the global features of the funnel-like surface change with the range of 
the potential, and to relate these changes to the dynamics of structural relaxation. We 
find that the landscape becomes rougher and less steep as the range of the potential 
decreases, and that relaxation paths to the global minimum become more complicated. 

1 Introduction 

Structural relaxation plays a key role in a diverse range of problems in chemical physics, 
including protein folding, glass formation, and the observation of "magic number" peaks 
in the mass spectrometry of rare gas clusters. The dynamic evolution of such systems is 
determined by the potential energy surface (PES) generated by the interactions between 
their constituent particles. Quite often one wants to find the structure and physical 
properties of a (macro) molecule or cluster, by which it is usually meant the properties of 
the global minimum on the PES, or, equivalently, the properties at zero Kelvin. However, 
the dynamics of a system at temperatures or energies above which it can escape from the 
global minimum depend on larger regions of the PES, the topology and topography of 
which determine the precise behavior. When considering the wider features of the PES in 
this way, it has become usual to refer to the PES as the "potential energy landscape" . 

One can also consider the free energy landscape, a temperature-dependent function 
which incorporates the entropy. For example, in protein folding such a landscape can be 
defined either as a function of the protein configuration by averaging the free energy over 
all solvent coordinates, or as a function of distance from the folded state in terms of a 
similarity parameter.^ 
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In recent years, much understanding has been gained in a number of fields by relating 
structural and dynamical properties to the underlying PES. For example, many years ago 
Levinthal pointed out the apparent contradiction between the astronomical number of 
possible configurations that a protein can adopt and the rapidity with which it finds the 
biologically active structure when it folds 12121 The "paradox" is resolved by realizing that 
efficient folding is only possible when the potential energy landscape is dominated by a 
funnel, i.e. consists largely of convergent kinetic pathways leading down in energy towards 
the required structured The precise features of a funnel may vary, but the native state 
must be thermodynamically stable at temperatures or energies where the dynamics are 
fast enough for the system to be able to explore the landscape and find itP The native 
state is destabilized if there are structurally distinct states of low energy which can act as 
kinetic trapsP Hence, a pronounced global minimum encourages efficient folding!^' 

The potential energy landscape also plays an important role in determining the be- 
havior of bulk liquids. Angell has proposed a widely used scheme in which liquids are 
classified from "strong" to "fragile" P A strong liquid is characterized by a viscosity whose 
temperature dependence follows an Arrhenius relationship (oc expL4/T]). These are often 
liquids with open network structures like water and Si02, whereas fragile liquids tend to 
have more isotropic interactions. AngelP^and Stillingei^ have described the general fea- 
tures of the energy landscapes that might be expected to characterize the two extremes. 
In a recent study, Sastry et al. have investigated the role of different regions of the land- 
scape in the process of glass formation in a model fragile liquid.^ They find that as the 
temperature of the liquid is decreased, the system samples regions with higher barriers, 
and on further cooling it samples deeper minima and non-exponential relaxation sets in. 

Another way that an energy landscape can be classified is as "sawtooth-like" or 
"staircase-like" depending on the energy difference between minima relative to the barriers 
which separate themf^E] For example, the "structure-seeking" properties of the (KC1)32 
cluster (i.e. its ability to find a rock salt structure even when cooled rapidly) can be at- 
tributed to downhill barriers which are low compared to the potential energy gradient 
towards crystalline minima, as in a staircase. 

In order to characterize an energy landscape, it is necessary to make a survey of its 
important features: minima, transition states and pathways. Since the number of such 
features increases at least exponentially with the number of particles in the system^ it is 
impractical and undesirable to catalogue them all for large systems. Consequently, existing 
studies have usually concentrated on analyzing what is hoped to be a representative sample 
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of minima and transition states U2llfil In this study we examine in detail the landscape of 
the 13-atom Morse cluster (M43), which is large enough to possess a complex PES, but is 
small enough for us to make a nearly exhaustive list of its minima and transition states. 
This model system is especially interesting because the energy landscape is dominated by a 
funnel, and the potential contains one parameter which allows us to adjust the complexity 
of the PES. Previous studies^^^ have shown that potential energy surfaces are simpler for 
short-ranged potentials, and the effects of the range on the morphology of global minima of 
atomic cluster j ^ 21 - 1 * and the stability of simple liquids-^^ have already received attention. 
The range of the potential also affects phase behavior: in a study of 7-atom Morse clusters, 
Mainz and Berry found that liquid-like and solid-like phase coexistence is less distinct when 
the range of attraction is longerP^ 

In this paper, we concentrate on finding useful ways to characterize and visualize a 
complex PES, and in the Summary we comment on how the range of the potential is likely 
to affect the relaxation properties of the cluster. We are currently using the data collected 
in this study to perform master equation dynamics on the system to address relaxation in 
detail. 

2 Exploring the Landscape 

The Morse potential can be written in the form 

V = Y, Vif, Vij = e" (1 - r ^ /rc) [ e '( 1 - r w/ r «> - 2]e, (1) 

i<j 

where is the distance between atoms i and j. e and r e are the dimer well depth 
and equilibrium bond length, and simply scale the PES without affecting its topology. 
They can conveniently be set to unity and used as the units of energy and distance. 
p is a dimensionless parameter which determines the range of the inter-particle forces, 
with low values corresponding to long range. Physically meaningful values range at least 
from p = 3.15 and 3.17 for sodium and potassium® to 13.62 for Cgo molecules When 
p = 6, the Morse potential has the same curvature as the Lennard-Jones potential at the 
minimum. 

The first step in characterizing the PES is to map out the local minima and the 
network of transition states^ and pathways that connects them. The eigenvector- following 
technique^H22l can efficiently locate transition states (first order saddles) by maximizing the 
energy along a specified direction, while simultaneously minimizing in all other directions. 
The minima connected to a given transition state are defined by the steepest descent paths 
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commencing parallel and antiparallel to the transition vector (the Hessian eigenvector with 
negative eigenvalue) at the transition state. Although eigenvector-following can also be 
used for these minimizations, the pathways are not necessarily the samej^ and may even 
lead to a different minimum. Since both the pathways and the connectivity are of interest 
here, we use a steepest descent technique for minimizations, employing analytic second 
derivatives, following Page and McIverP^ 

Our algorithm for exploring the PES is similar to that used by Tsai and Jordan in a 
study of small Lennard-Jones and water clusters.^ Starting from a known minimum: 

1. Search for a transition state along the eigenvector with the lowest eigenvalue. 

2. Deduce the path through this transition state and the minima connected to it. 

3. Repeat from step^beginning antiparallel to the eigenvector, and then in both direc- 
tions along eigenvectors with successively higher eigenvalues until a specified number, 
n ev , of directions have been searched uphill. 

4. Repeat from step^until n ev modes of all known minima have been searched. 

By taking steps directly between minima, this method avoids wasting time on intra-well 
dynamics. Other methods for exploring energy landscapes, such as molecular dynamics, 
can become trapped in local minima, especially at low temperature, where there is a wide 
separation in time scale between inter-well and intra-well motion. The chosen value of n ev 
clearly affects the thoroughness of the survey, although even if all (3N — 6) vibrational 
modes of an iV-atom cluster are searched, there is no guarantee of finding every minimum 
and transition state. In practice, the required computer time and storage demand that 
n ev be reduced for large p, since the complexity of the PES increases dramatically as the 
range of the potential decreases. However, one finds that searches from low- lying minima 
are more likely to converge in a reasonable number of iterations, so the above algorithm 
was augmented with searches along further eigenvectors of lower-energy minima. We are 
confident that the databases generated for p = 4 and 6 are nearly exhaustive, and although 
those for higher values of p are necessarily less complete, this approach still allows us to 
map out the PES fairly comprehensively. 

Details of the searches and the resulting databases for p = 4, 6, 10 and 14 are summa- 
rized in Tabled The dramatic rise in the number of minima and transition states found 
as the range of the potential decreases is the first indication of the increasing complexity 
of the PES.^ The remainder of this paper investigates in more detail the nature of these 
changes and some useful ways of characterizing the landscapes. 
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3 Topological Mapping 



When trying to describe an energy "landscape" , one has already been forced to use termi- 
nology appropriate to a surface in three-dimensional space, and pictorial representations 
are usually restricted even further to two dimensions. Visualizing a 3iV-dimensional object 
directly in such a way has obvious limitations, yet it is appealing to have an idea of "what 
the surface looks like". 

One helpful way of doing this is to use topological mapping to construct a disconnec- 
tivity graph, as applied to a polypeptide by Becker and KarplusP^The analysis begins by 
mapping every point in configuration space onto the local minimum reached by following 
the steepest descent pathf^ Thus, configuration space is represented by the discrete set 
of minima, each of which has an associated "well" of points which map onto it. Although 
this approach discards information about the volume of phase space associated with each 
minimum, the density of minima can provide a qualitative impression of the volumes 
associated with the various regions of the landscape. 

At a given total energy, E, the minima can be grouped into disjoint sets, called basins 
("super basins" in Becker and Karplus' nomenclature), whose members are mutually ac- 
cessible at that energy. In other words, each pair of minima in a basin are connected 
directly or through other minima by a path whose energy never exceeds E, but would 
require more energy to reach a minimum in another basin. At low energy there is just one 
basin — that containing the global minimum. At successively higher energies, more basins 
come into play as new minima are reached. At still higher energies, the basins coalesce as 
higher barriers are overcome, until finally there is just one basin containing all the minima 
(provided there are no infinite barriers). 

The disconnectivity graph is constructed by performing the basin analysis at a series 
of energies, plotted on a vertical scale. At each energy, a basin is represented by a node, 
with lines joining nodes in one level to their daughter nodes in the level below. The choice 
of the energy levels is important; too wide a spacing and no topological information is left, 
whilst too close a spacing produces a vertex for every transition state and hides the longer 
range structure of the landscape. The horizontal position of the nodes is arbitrary, and 
can be chosen for clarity. In the resulting graph, all branches terminate at local minima, 
while all minima connected directly or indirectly to a node are mutually accessible at the 
corresponding energy. 

The disconnectivity graphs for M13 with p = 4 and 6 are plotted on the same scale 
in figure ^ We have chosen a linear energy spacing of one well depth, which is an effec- 
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tive compromise between the points raised above. Both trees are typical of a funnel-like 
landscape: as the energy is lowered, minima are cut off a few at a time with no secondary 
funnels, which would appear as side branches. A large upward shift in the energy range 
of the minima is apparent on increasing p from 4 to 6, due to the increase in the energetic 
penalty for strain and a decrease in the energetic contribution from next-nearest neighbors 
as the range of the potential decreases.^ An increase in barrier heights is also revealed by 
the somewhat longer branches at p = 6. Because of the large number of minima involved 
in the databases for p = 10 and 14, the disconnectivity graphs are too dense to illustrate, 
but we shall see in the numerical analysis of the next section how the trends develop. 

The concepts involved in the disconnectivity graph have much in common with the 
"energy lid" description of Sibani et alP^ in which minima are grouped together if they 
are connected by paths never exceeding a particular energy (the "lid"). These authors 
plotted a tree with a time axis, on which nodes represent the time when groups of minima 
first come into equilibrium. 

The term "basin" has been used with a somewhat different meaning by Berry and 
coworkers! 1 ^ 34 * In this definition, a basin consists of all minima connected to the basin 
bottom by a monotonic sequence, i.e. a sequence of connected minima with monotonically 
decreasing energy. This definition contrasts with that of Becker and Karplusp^l because it 
is independent of the energy, and actually has a lot in common with the notion of a funnel. 
Although the word "funnel" may conjure up a misleading image when the surface is rough 
or shallow in slope, we will use it in this context to avoid confusion with the previous 
definition of a basin as a set of mutually accessible minima at a given energy. The funnel 
terminating at the global minimum is denoted the primary funnel, whilst adjoining side 
funnels are termed secondary. It should be noted that this definition permits a minimum 
to belong to more than one funnel via different transition states. The significance of 
dividing the landscape in this way is that inter-funnel motion is likely to occur on a slower 
time scale than inter-well flow, 1^231 so funnels constitute the next level in a hierarchy of 
landscape structure. Sufficiently deep or voluminous secondary wells can act as trapsP^ 
A striking example is the cluster of 38 Lennard- Jones atoms, whose truncated octahedral 
global minimum was only found quite recentl ^9 l 35 l because of the much larger secondary 
funnel associated with a low-lying icosahedral structured 

As the first line of Table [21 shows, for p = 4 the landscape of M13 is a perfect funnel: 
all minima lie on monotonic sequences terminating at the global minimum. At higher 
values of p a small fraction of minima lie outside the primary funnel, and although they 
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technically constitute secondary funnels, they represent a very small proportion of the 
phase space. We will now see how the characteristics of the primary funnel evolve as the 
range of the potential is decreased. 

4 Properties of the Landscape 

The remainder of Table El lists some global properties of the landscape at four values 
of p. Some of the trends are straightforward to understand. For example, defining \>i 
as the geometric mean of the normal mode frequencies at minimum i, the average of 
this quantity over the database of minima, {P) m , rises monotonically with p because of 
the increasing stiffness of shorter-ranged potentials. The average of the transition state 
imaginary frequency, z^ m , increases less rapidly in magnitude, and levels off at high p, 
indicating that the transition regions are flatter relative to the well bottoms than at low 
P- 

The increasing energy, E gm , of the global minimum was noted in the previous section, 
and the table shows that this increase is accompanied by a decreasing gap AE gSuP to the 
second lowest minimum. The striking drop in AE gap when p reaches 14 is due to a change 
in morphology of the second lowest structure, as illustrated in figure El To see why this 
happens, it is helpful to decompose the potential energy into the following contributions!^ 

V = — n nn + Entrain + -^nnn, (2) 

where n nn is the number of nearest neighbor interactions, i.e. the number of pairs lying 
closer than a value ro (taken here to be 1.15r e ), and the strain energy and non-nearest 
neighbor contributions are defined by 

Strain = £ [Vy + 1], (3) 

nj<r 

£nnn= £ V ly (4) 

i<j 

n nn and -Entrain are more sensitive properties of the structure than E nnn , and so the lowest 
energy cluster is determined by a balance between maximizing n nn and minimizing -Entrain- 
The icosahedron [figure |3f a)] is the global minimum for all four values of p considered here 
because it has the largest number of nearest neighbors (n nn = 42). However, the large 
value of n nn is at the expense of considerable strain. As E stra [ n is the energetic penalty for 
nearest-neighbor distances deviating from r e , it increases rapidly for strained structures 
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as the pair-potential well narrows at larger p. E nnn is also sensitive to p; it decreases as 
the range of the potential decreases. 

The upward trends in figure |2] are caused by the changes in -Entrain and E nnn . For 
p < 13.90 the second lowest minimum is a defective icosahedron in which one vertex 
has been removed and one face is capped [figure EJb)]. The removal of a vertex allows 
the strain in the icosahedron to relax, and so the energy rises less steeply than for the 
icosahedron and A-Eg ap falls. However, decahedral clusters are intrinsically less strained 
than icosahedral ones, and at p = 13.90 the decahedron [figure E(c)] , which for lower 
p is a transition state, becomes the second lowest minimum. In fact, for p > 14.77 
the decahedron is the global minimum, although this value of p may be too large to be 
observed in chemical systems. The change in the order of the stationary point arises from 
a delicate balance between -Entrain an d -Ennn- The vibrational mode of the decahedron 
with the lowest Hessian eigenvalue is a twist about the C5 axis. This motion strains the 
structure, but brings non-nearest neighbors closer. At high p, the increased strain wins, 
causing the energy to rise and giving a minimum, whereas for longer-ranged interactions 
the non-nearest neighbors lower the energy, giving a saddle. 

The decreasing A-E gap indicates a local flattening of the PES at higher p. This effect 
extends beyond the vicinity of the global minimum to the whole landscape, as can be seen 
from the energy distributions of minima shown in figure As p increases, the energy 
distribution shifts upwards and becomes narrower, and for p = 10 and 14 it develops 
two sharp peaks at —33 and —34. At high values of p, E nnn becomes small, and the 
energetic penalty for strain is large. Decomposition of the energy according to equation 
121 reveals that the peaks in the distributions correspond to low-strain structures with 33 
and 34 nearest neighbors. Low strain can arise from two structural motifs: close packing 
or polytetrahedral packing (without pentagonal bipyramids). It is not easy to classify 
such a small cluster according to these schemes, but it is worth noting that the radial 
distribution function, taken over all the minima, develops a \[2 signature as p increases, 
which is characteristic of close packing.^ 

A more quantitative measure of the slope of the PES is provided by the energy differ- 
ence between pairs of connected minima, AEf on (where i labels the connecting transition 
state, or, equivalently the pathway). As Table EJ shows, the average of this quantity over 
the pathways drops off quickly as p increases from 4. AEf on is the difference between 
the uphill and downhill barriers 6" p and bf own defined by transition state i and the two 
minima it connects. Although the average over the pathways of the uphill barrier, (6 up ) p , 
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decreases as the range of the potential decreases, (6 down ) p increases, i.e. the barriers that 
must be overcome for structural relaxation towards the global minimum are larger; the 
flattening of the funnel is accompanied by roughening. 

Given the dramatic increase in the number of stationary points as the range of the 
potential decreases, and that the volume of accessible phase space will be reduced as the 
long range attraction is squeezed out, we must expect some change in the nature of the 
individual pathways between minima and their organization on the landscape. Defining Di 
as the separation in configuration space of the two minima connected by transition state 
i, Table [21 shows, as we might expect, that connected minima are on average closer when 
the potential is short-ranged. This effect is accompanied by a decrease in the average 
of the integrated path length, Si. It is interesting to see how the individual pathways 
are organized into routes to the global minimum. We have calculated the shortest path 
from each minimum to the global minimum, as measured by the total integrated path 
length Sf m (the path with fewest steps between minima is generally longer). The average 
of Sf m is fairly insensitive to p, whilst the average of the number of steps along the 
corresponding pathways, nf m , increases. Thus, on average, the path for relaxation to the 
global minimum does not increase significantly in length, but becomes more rugged as 
more transition states must be crossed. Whereas every minimum at p = 4 can reach the 
global minimum in either 1 or 2 steps, as many as 5 may be required at p = 14. Table 01 
shows how the minima are distributed over n? m , giving some insight into the connectivity 
of the landscape. The number of minima with nf m = 1 tells us how many transition states 
are connected directly to the global minimum. The values are remarkably high, especially 
as permutational isomers are not included. Interestingly, the number of minima does not 
increase continuously as the sequences branch out from the global minimum (as one might 
expect in a funnel), but tails off quite gently. 

An intuitive explanation for the constancy of {S gm ) and the increase in (n gm ) might 
be that paths are split into a larger number of sub-rearrangements. The number of atoms 
contributing to rearrangement i can be measured by the cooperativity index iV, = N/ji, 
where ji is the moment ratio of displacement, which is defined by^ 

(E£ Ms) - r a (i)P) 

where r a is the Cartesian position vector of atom a, and s and t denote the final and initial 
configurations in rearrangement pathway i. Table [21 shows that the average value of iVj is 
almost independent of p. In fact the distribution of Ni (from 1 to N) is remarkably similar 
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for all four databases. This result contrasts with statistics previously obtained for the 
larger clusters LJ55 and (C6o)s5, which showed that cooperative (high iVj) rearrangements 
are less likely for (C6o)55> where the range of the potential corresponds to p « 14)221 It 
is possible that a 13-atom cluster is too small to support localized sub-rearrangements in 
this way. 

5 Summary 

We have performed a comprehensive survey of the potential energy landscapes of the 13- 
atom Morse cluster for four values of the range parameter using systematic eigenvector- 
following searches. The landscapes were then characterized in detail using disconnectivity 
graphs, funnel analysis, and a selection of parameters that provide insight into the topology 
and topography. We have described and rationalized the changes in the landscape as the 
range of the potential is varied over a physically meaningful range. 

The trends displayed in Table Hand the above discussion are underlined by the plots of 
representative monotonic sequences in figure El The overall classification of the potential 
energy landscape is that of a funnel, but one which becomes flatter and rougher as the 
range of the potential decreases. This change is accompanied by a general increase in 
complexity of the surface in terms of the number of minima and transition states and in 
the number of steps required to reach one minimum from another. 

Previous studies of model potential landscapes^ have shown that relaxation from high- 
energy configurations to the global minimum is most efficient when the PES has a large 
potential energy gradient towards the global minimum with low downhill barriers, and 
lacks secondary funnels which act as kinetic traps. On this basis we would expect M13 
to relax most easily when the the range of the potential is long, in spite of the fact that 
the frequency of intra-well vibrational oscillations decreases as the potential becomes less 
"stiff" at fixed values of e and r e (see Table |2J). Low values of the range parameter p are 
therefore likely to produce "structure-seekers" whereas high values will tend to produce 
"glass- formers" , reflecting a continuous change from a staircase-like to a sawtooth-like 
landscape. 

Equipped with an understanding of the potential energy landscape and its dependence 
on the range of the potential, we have applied the master equation approach to investigate 
the dynamics of structural relaxation in M13 . This work enables us to probe in detail the 
flow of probability between individual minima in an ensemble of clusters as they approach 
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the equilibrium distribution, and the results will be described in a separate publication. 

Acknowledgments 

M.A.M. and D. J.W. gratefully acknowledge financial support from the Engineering and 
Physical Sciences Research Council, and the Royal Society, respectively. M.A.M. would 
also like to thank Gonville & Caius College, Cambridge for a grant towards attending the 
workshop on energy landscapes at the Telluride Summer Research Center in July 1997. 

References 

[1] J. D. Bryngelson, J. N. Onuchic, N. D. Socci and P. G. Wolynes, Proteins: Struct., 
Func. and Gen. 21, 167 (1995). 

[2] C. Levinthal in Mossbauer spectroscopy in biological systems, proceedings of a meeting 
held at Allerton House, Monticello, Illinois, edited by P. DeBrunner, J. Tsibris and 
E. Munck, p. 22. University of Illinois Press (1969). 

[3] Ref. |21 is not readily available in libraries, but it can now be downloaded from URL 
http : //brian . ch . cam . ac . uk/~mark/levinthal/levinthal . html 

[4] P. E. Leopold, M. Montal and J. N. Onuchic, Proc. Natl. Acad. Sci. USA 89, 8721 
(1992). 

[5] P. G. Wolynes, Proc. Natl. Acad. Sci. USA 93, 14249 (1996). 

[6] J. P. K. Doye and D. J. Wales, J. Chem. Phys. 105, 8428 (1996). 

[7] A. Sali, E. Shakhnovich and M. Karplus, J. Mol. Biol. 235, 1614 (1994). 

[8] P. Amara and J. E. Straub, J. Phys. Chem. 99, 14840 (1995). 

[9] C. A. Angell, J. Non-Cryst. Solids 131-133, 13 (1991). 

[10] C. A. Angell, Science 267, 1924 (1995). 

[11] F. H. Stillinger, Science 267, 1935 (1995). 

[12] S. Sastry, P. G. Debenedetti and F. H. Stillinger, Nature 393, 554 (1998). 

[13] R. E. Kunz and R. S. Berry, J. Chem. Phys. 103, 1904 (1995). 



11 



[14] K. D. Ball, R. S. Berry, R. E. Kunz, F.-Y. Li, A. Proykova and D. J. Wales, Science 
271, 963 (1996). 

[15] F. H. Stillinger and T. A. Weber, Phys. Rev. A 25, 978 (1982). 

[16] J. P. K. Doye and D. J. Wales, J. Phys. B 29, 4859 (1996). 

[17] M. R. Hoare and J. Mclnnes, Faraday Discuss. 61, 12 (1976). 

[18] P. A. Braier, R. S. Berry and D. J. Wales, J. Chem. Phys. 93, 8745 (1990). 

[19] J. P. K. Doye, D. J. Wales and R. S. Berry, J. Chem. Phys. 103, 4234 (1995). 

[20] J. P. K. Doye and D. J. Wales, J. Chem. Soc, Faraday Trans. 93, 4233 (1997). 

[21] J. P. K. Doye and D. J. Wales, Science 271, 484 (1996). 

[22] D. T. Mainz and R. S. Berry, Mol. Phys. 88, 709 (1996). 

[23] P. M. Morse, Phys. Rev. 34, 57 (1929). 

[24] L. A. Girifalco and V. G. Weizer, Phys. Rev. 114, 687 (1959). 

[25] L. A. Girifalco, J. Phys. Chem. 96, 858 (1992). 

[26] J. N. Murrell and K. J. Laidler, J. Chem. Soc, Faraday Trans. 64, 371 (1968). 

[27] J. Pancir, Coll. Czech. Chem. Comm. 40, 1112 (1974). 

[28] C. J. Cerjan and W. H. Miller, J. Chem. Phys. 75, 2800 (1981). 

[29] D. J. Wales, J. Chem. Phys. 101, 3750 (1994). 

[30] M. Page and J. W. Mclver, J. Chem. Phys. 88, 922 (1988). 

[31] C. J. Tsai and K. D. Jordan, J. Phys. Chem. 97, 11227 (1993). 

[32] O. M. Becker and M. Karplus, J. Chem. Phys. 106, 1495 (1997). 

[33] P. Sibani, J. C. Schon, P. Salamon and J. Andersson, Europhys. Lett. 22, 479 (1993). 

[34] R. S. Berry and R. Breitengraser-Kunz, Phys. Rev. Lett. 74, 3951 (1995). 

[35] J. Pillardy and L. Piela, J. Phys. Chem. 99, 11805 (1995). 

[36] D. J. Wales, M. A. Miller and T. R. Walsh, Nature (in press 1998). 

12 



[37] J. P. K. Doye and D. J. Wales, Chem. Phys. Lett. 247, 339 (1995). 

[38] N. W. Ashcroft and N. D. Mermin, Solid State Physics, Holt, Rinehart and Winston, 
New York (1976). 

[39] F. H. Stillinger and T. A. Weber, Phys. Rev. A 28, 2408 (1983). 



13 



Table 1: Details of the databases for M13 at four values of the range parameter p. n cv is 
the minimum number of eigenvectors of each minimum searched for a transition state, and 
n s is the average number of searches from each minimum. n m i n and n ts are the numbers 
of minima and transition states found. 



p 


4 


6 


10 


14 




15 


6 


3 


2 


n s 


31.3 


13.0 


7.0 


7.5 




159 


1439 


9 306 


12 760 


rhs 


685 


8 376 


37499 


54439 
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Table 2: Some properties of the potential energy landscape of M13 at four values of the 
range parameter p. All dimensioned quantities are tabulated in reduced units. n m \ n is 
the number of minima, of which n p f lie in the primary funnel. E gm is the energy of the 
global minimum, with the next-lowest energy structure lying AE ga _ p higher, is the 
geometric mean normal mode frequency at minimum i and u™ 1 is the imaginary frequency 
at transition state i. 6" p is the larger (uphill) barrier height between the two minima 
connected by transition state i, bf own is the smaller (downhill) barrier, and AEf on is the 
energy difference between the minima, so that bf = bf own + AEf on . Si is the integrated 
path length between the two minima connected by transition state i, Di is their separation 
in configuration space, and iVj is the cooperativity index of the rearrangement (defined in 
the text). nf m is the smallest number of steps from minimum i to the global minimum, and 
Sf m is the integrated length of this path. (• • -) m , (• • -) ts and (• • -) p indicate averages where 
the index runs over minima, transition states, and non-degenerate pathways (i.e. pathways 
not merely connecting permutational isomers) respectively. 
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Table 3: The distribution of the number of steps n gm lying on the shortest path from local 
minima to the global minimum at four values of the range parameter p. 
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Figure 2: Correlation diagram for some low- lying structures (see ngureEJ): the icosahedron 
(Ih), the decahedron (-D5/J and the lowest-energy defective icosahedron (C s ). Dashed lines 
indicate regions where the structure is not a minimum: D$h becomes a transition state 
and C H becomes a second order saddle. 
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(a) (b) (c) 

Figure 3: Structures discussed in the text: (a) the icosahedron (i/j), (b) the lowest-energy 
defective icosahedron (C s ), and (c) the decahedron (D^)- 
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Figure 4: Energy distribution of the minima for four values of the range parameter p. In 
each case, the energy of the global minimum is indicated by an arrow. 
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Figure 5: Example monotonic sequences leading to the global minimum for three values of 
the range parameter p. S is the integrated distance along the reaction path from the global 
minimum. Minima are indicated by filled circles, and transition states by open circles. The 
plots demonstrate a number of features discussed in the text: the general increase in energy 
of the minima, the decreasing gap to the global minimum, the increasing barrier heights, 
the shorter rearrangements, and the decreasing gradient towards the global minimum as 
the range of the potential decreases. 
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